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Abstract 

Recent data on the Gottfried sum and less recent ones on the pion-nucleon 
sigma term seem to disagree with naive parton-based expectations on the light 
(i.e., up, down and strange) quark content of the nucleon. We show that these 
discrepancies are resolved if nonperturbative contributions are included in the 
analysis of the data. These appear both in the computation of matrix elements 
of operators, and in their QCD scale dependence, and depend strongly on the 
quantum numbers of the given state. 
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1 Introduction 



In recent times a set of high-precision measurements of the nucleon structure func- 
tions, both in the polarized and unpolarized case, has lead to results in striking 
disagreement with naive parton model expectations. The first of these to come [|l| 
has spawned a vast theoretical literaturep|, the so-called nucleon spin crisis; it is now 
clear that the experimental result can be fully explained and understood only in 
a nonperturbative framework, and, more specifically, that it may be a first hint of a 
nonperturbative effect which survives the high-energy limit and is thus visible in the 
deep-inelastic regime. 

More recent results on the Gottfried sum rule [^[Q have led to question even 
the standard assumptions on the up and down (unpolarized) flavor content of 
the nucleon. Combined with the long-standing sigma term puzzle 0, which hints to 
an unusual strange quark content, these results seem to suggest that the light-flavor 
structure of the nucleon is rather different from what the naive intuition based on a 
combination of the non-relativistic quark model with the naive parton model would 
suggest. 

It seems thus that naive expectations have to be reconsidered. These are typically 
expressed as Ansatze on the form of the quark distributions qf{x), which in the 
parton model provide the probability of a quark of flavor i to carry a fraction x of 
the momentum of the baryon k, and are related to the nucleon structure functions. 
The quark distributions are split into a valence and part: 

q\x) = q'^{x) + q'^{x); q\x)=ql{x). (1) 

In QCD this decomposition has only meaning as a definition: the quark distributions 
are related to the expectation value of certain operators in such a way that the 
structure function F2{x) is related to them by[0 

F^{x) = xY.el(qK^) + q^{x)), (2) 

i 

where the sum runs over all the quark flavors and Cj are electric or weak charges 
according to whether the target is probed with electrons or neutrinos. The quark 
and antiquark components for various flavors can be extracted by taking suitable 
linear combinations of structure functions for neutrino, antineutrino and electron 
scattering off different targets. QCD then leads to predict the scale dependence of 
these quantities. 

However, it is usually assumed that the quark distributions can be actually viewed 
as expressing the quark content of the given hadron, and that this picture provides 
information on the nucleon structure in that it suggests the form of the distribution 
themselves. This assumption is based on the expectation that the matrix elements 
of operators which provide the quark-parton distributions (and which can only be 
written down in the light-cone formalism 0) can be identified, at least up to some 
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unknown renormalization, with the quark operators of the naive quark modeL Such 
an assumption is totally unjustified within QCD; it is however quite successful phe- 
nomenologically, and it is customarily used in analyzing the experimental data. 

A simple example of application of this wisdom is the computation of the ratio 
R{x) = as x — > 1. The naive argument goes as follows: as x — > 1 only 

the valence quarks survive, because if a quark is carrying the whole of the nucleon's 
momentum it must be a valence one. Furthermore, if one single quark is carrying the 
quantum numbers of a nucleon, by isospin, it ought to be a -u quark for the proton 
and a d quark for the neutron. Thus lim^.^! = e^/e^ = 1/4. This prediction 
is spectacularly borne out by the data 0],|]TO|, while it is clearly based partly on 
symmetry (isospin) and partly on naive guesswork. 

More in general, the quark distributions (|I|) are usually taken to satisfy the fol- 
lowing assumptions: 

1. The sea and valence components separately satisfy isospin symmetry: 

ui = d:, u: = di- ui = d:, u: = di. (3) 

2. The sea component in each baryon is SU(2) symmetric: 

ui = di- u: = d:. (4) 

3. The content of Zweig suppressed flavors, such as strange quarks in the nucleon, 
is small. This implies relations like 

{p\ss\p) ~ {n\ss\n) ^ 0. (5) 

It turns out that the recent data on the Gottfried sum rule put either (or both) as- 
sumptions 1 and 2 into question, and the long-standing sigma term puzzle challenges 
assumption 3. In the sequel, we will discuss both these points in turn. We shall show 
that the naive picture can actually be recovered, provided nonperturbative QCD ef- 
fects are taken into account, both in the QCD evolution and in the determination of 
the symmetry structure of matrix elements of operators. 



2 The puzzle of the Gottfried sum 

The question of the validity of assumptions 1 and 2 can be condensed into an analysis 
of the Gottfried sum So, defined by 

s,./,,£M^^. (6) 
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According to Eq.(|^), the first moment of F2 corresponds to tlie sum of tlie electric 
cfiarges squared of all partons present in the target, 



i i 

where we have defined 

ni= dx {qi{x) + qi{x)) . (8) 



This is expected to be a divergent quantity for a single nucleon, due to the fact 
that at small x the nucleon's momentum is subdivided among an infinite number of 
quarks; however this divergence (i.e., the small-x behavior of all q\x)) is expected 
to be universal (on the basis of Regge theory) [p!| , essentially due to the flavor- 
independence of the strong interaction, hence So Eq-(§D is expected to be finite. 

The symmetry assumptions (1),(2) amount to the further conjecture that not only 
in the a; — > limit, but actually for all x the sea is universal, hence only valence quarks 
can contribute to nonsinglet combinations such as Eq.(^. Taking only the valence 
component into account one arrives thus at the so-called Gottfried sum ruleJTT| ]T^: 



4 14 1 

Sg = g<,„ + - g<,„ " g^^ (S) 

j - 1) - <J = I- (10) 

However Sq has now been measured^ as 

= 0.24 ± 0.016 at 4 GeV^. (11) 

Even though the extraction of this value involves certain assumptions, e.g. an extrap- 
olation of the measured structure function to the full range of x, it seems exceedingly 
unlikely that a value compatible with the naive expectation Eq.(p^) could be ob- 
tained. As a matter of fact, an improved data analysis ^ has given the new value 

5g = 0.258 ±0.017 at 4 GeV^. (12) 

but has also shown evidence for shadowing in deuterium. If shadowing effects (which 
cannot be measured precisely as yet) were included, the value of Sg as currently 
determined (Eq.(|lT]) or (0)) should be brought down by an amount which is the- 
oretically estimated |[T3[] to be ASg ~ 0.2 — 0.4. The value Eq. (|Tl|) seems thus a 
conservative estimate of the value of Sq, which might turn out to be even lower. 

It is thus necessary to reconsider assumptions 1 and 2. If we make no assumption, 
the Gottfried sum is generally given by 
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^{riu + rid) + ^(^« - nd) 



(13) 

1=1 
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where the subscript indicates that the isotriplet part must be taken, i.e. the difference 
of the expression in square brackets between a proton and a neutron. Assumption 1 
(isospin symmetry) imphes [nu+nfi\i=i = 0, whereas assumption 2 (seaffavor symme- 
try) imphes that [nu — nd]i=i be identified with its valence value, i.e., [nu — nd]i=i = 1. 
The data could be explained either by assuming that the neutron sea is bigger than 
the proton one (thus violating isospin), or that the proton sea and separately the 
neutron sea are not flavor neutral, but rather they have a flavor asymmetry that 
anticorrelates to that of the valence component (thus violating sea flavor symmetry). 

In order to fix the relative amount of violation of these two symmetries one needs 
an independent measurement of a different linear combinations of the quantities en- 
tering Eq. ([l3|) . It can be shown [p^ that the required information is provided by 



the isotriplet part ac of the so-called nucleon cr-term, defined as 

ac = o-/=i = {N\muuu + ■mddd\N)j=i, (14) 



where rrii are current quark masses. This is because ac satisfies [jT^ [15 



"^G = ^ [{mu + md){nu + Ud) + (m„ - nid){n^ - nd)]^^^ = Mp - (15) 

which is the analogue of the Gottfried sum Sq, with square charges replaced by 
masses. Using this extra piece of information one gets to the conclusion that 

l[nu + nd]i=i ^ ±0.03 (16) 

^[n^-nd]i=i-l ^ -0.35, ±0.1 (17) 

i.e., sea flavor neutrality is broken considerably more strongly than isospin. 

Just as in the well-know case of the "spin puzzle" [00 it turns out that the 
experimental result Eq.([T^ ) - ([l7|) can be easily accommodated in effective models of 
the nucleon, such as the Skyrme or the bag model |]T1|. Thus, again like in the spin 



case, one may ask whether the data are not simply telling us that the "naive parton 
wisdom" is simply not applicable in this channel. 

To make this more precise, one should first however investigate what QCD does 
tell us. Any nonsinglet combination of first moments of structure functions F2{x), 
such as Sg Eq-(§|) acquires a two-loop scale dependence via the Altarelli-Parisi 
equationspl. Indeed, it is clear that at one loop the quark number does not evolve, 
because at one loop the only elementary process (in the nonsinglet channel) is gluon 
radiation by a quark line (Fig. 1) which clearly leaves the total number of quarks 
(albeit not their momentum distribution) unchanged. At two loops, however, the 
emitted gluon can in turn produce a quark- ant iquark pair through the diagram of 
Fig. 2, thereby increasing the total quark plus antiquark number measured by the first 
moment of F2{x). At first sight it would seem that this again should not contribute 
to the evolution of nonsinglet quantities, because the probability for pair emission 
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is independent of the emitted flavor. However, this is not entirely correct because 
if the emitted pair has the same flavor as the original quark the final state must be 
antisymmetrized with respect to the two identical quarks, while this is of course not 



the case if the flavor of the radiated pair is different 1 16]. This is enough to gener- 



ate a tiny difference between the splitting function "P^ for flavor- diagonal emission 



(i.e., the case when the emitted pair has the same flavor as the original quark) and 
the splitting function Vqq^ for flavor- nondiagonal emission (i.e., the case when the 
emitted pair has a different flavor). 



It follows that Sg acquires a scale dependence |[16|1; this is given by the two-loop 
Altarelli-Parisi equation]^ 

^( 

dr 

where ® indicated convolution of the splitting function and the quark distributions, 
NS indicates any nonsinglet combination of quark distributions, and the splitting 
function responsible for the nonsinglet evolution is given by 

Q = V°- V^^ (19) 

in terms of the diagonal and nondiagonal splitting functions. It follows that the 
evolution of Sq Eq-(i) is given by 

^^SG{t) = {Q\q + Q\q)SG{t), (20) 

in terms of the first moments Qgg, Q\q of the various sphtting functions. The physical 
interpretation of Eq. ([T8| ) (pO| ) (which is derived e.g. in Ref.|^) is straightforward: the 
nonsinglet evolution is due to the difference in probability for flavor- diagonal and 
flavor-nondiagonal pair production; if the latter is greater than the former then the 
nonsinglet quark number measured by Sg decreases because the sea generated from 
the process of Fig. 2 tends to compensate the starting flavor asymmetry. 

The same results can be obtained via operator-product expansion methods. 



which are however much more cumbersome because only even moments of the (charge- 
conjugation even) combination of structure functions which appears in Eq.(^ have 
a direct interpretations in terms of matrix elements of twist-two local operators, 
whereas the odd moments (such as the first moment which defines Sg) have to be 
constructed by analytic continuation. Anyway, Sg is thus found to evolve perturba- 
tively according to 



Saif-i'), (21) 



where explicit perturbative computation up to two loops leads to 'j^'^^ /b = 0.01 [|T| 
It follows that perturbative QCD evolution, even though present, is entirely negli- 
gible from a quantitative point of view because if < 1 (as it must in order for 
perturbation theory to make sense) then the total variation of Sg is less than 1 %. 
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It seems thus that, just hke in the spin case, even if we assume that the "parton 
wisdom" should apply at some low scale, while the experimental result is taken at 
large Q^, perturbative QCD evolution cannot explain the discrepancy between the 
experimental value Eg . ([TT|) and the expectation Eq.(^). More specifically, even if 
we are willing to give up the naive expectation Eq.(|T^) we must accept the rather 
unpalatable conclusion that even at very low scales the quark-parton content of a 
nucleon is very different from its constituent quark content. 

There is, however, a way out of this impasse, which consists of allowing for the 
possibility that the perturbative evolution Eq.(pl|) be replaced by some nonpertur- 
bative evolution at low enough scales. Due to a nonperturbative mechanism this 
evolution could then be much larger than the perturbative one. Symmetry consider- 
ations suggest what the origin of such a mechanism may be |T^. If the interaction 
which governs the generation of the quark-antiquark sea is U(Nj-)-symmetric, then 
by Wigner-Eckhart's theorem the sea will be flavor symmetric. In such case there is 
no QCD evolution at all, and indeed Eq.([T8|) shows that all non-singlet scale depen- 
dence can be traced to a difference in diagonal and non-diagonal splitting functions 
which signals the violation of exact U(Nj) symmetry. The evolution is thus, crudely 
speaking, a measure of the amount of U(Nj) symmetry breaking. Now, it is well 
known that in QCD the axial (Nj) flavor symmetry is dynamically broken down 
to a SUa (Nj) by the axial anomaly. This is manifested in the mass spectrum of 
pseudoscalar mesons, since the singlet meson (the 1]') has a much larger mass than 
the octet mesons (pions and kaons) which can be viewed as Goldstone bosons. If such 
effect is incorporated in the Altarelli-Parisi equations it could thus potentially lead to 
large evolution. Notice that this can only be done in an effective way, since there is 
no known way of obtaining this pattern of symmetry breaking through purely pertur- 
bative computations (rather, nonperturbative effects such as instantons are thought 
to be at the origin of it). 

We can thus try to effectively include nonperturbative symmetry-breaking ef- 
fects by considering contributions to the evolution equations where a meson couples 
directly to the quark-partons. These can be viewed pictorially as being given by 
diagrams with the same topology as those of Fig. 2, but where the original quark and 
the emitted antiquark exchange many gluons in such a way that they can be viewed 
as a quark-antiquark bound states, thus leading to the diagram depicted in Fig. 3. 
It is clear that, qualitatively, this should give a contribution to the Altarelli-Parisi 
evolution which may potentially generate signiflcant evolution with the right sign as 
required to explain the observed effect Eq. (PD . Indeed, due to the larger mass of 
the 7]', the singlet meson emission should be energetically suppressed in comparison 
to octet meson emission. Because flavor-nondiagonal emission is entirely due to the 
latter process, whereas flavor- diagonal emission proceeds in part through the former, 
one expects that the diagonal process will be disfavored. If Eq.(|18D may be gen- 
eralized to meson emission, this should result in a negative value of its r.h.s., thus 
leading to a partial screening of the valence flavor asymmetry. Notice that the pro- 
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cess given by the elementary diagram of Fig.3 is now effectively a one-loop one, thus 
the corresponding evolution is potentially strong. 

It seems thus reasonable to investigate this option in a quantitative way [^]. 
The required generalization of the Altarelli-Parisi evolution proceeds in two steps: 
first, we must generalize the set of evolution equations so as to include the coupling 
of quarks to bound states, and then we must compute the new splitting functions 
which enter the generalized equations. Because we are only interested in nonsinglet 
evolution, we may concentrate on the pseudoscalar meson sector, since for all other 
mesons or baryons the flavor symmetry will be approximately exact, thus, by the 
above argument, no contributions to the evolution should arise. 

The main difference between the usual Altarelli-Parisi equations and the general- 
ized ones is then the possibility of flavor-nondiagonal transitions already at the one 
loop level, such as that given by the diagram of Fig. 3 when one starts with, say, a u 
quark and ends with a d quark and a vr"*". It is then easy to see that the nonsinglet 
evolution equations generalize (considering for simplicity flavor SU(2) only) as 

= Qqq(S)q+ Qqq(^q + Vq^(g)n, (22) 
J^q = Qqq®q+Qqq®q-Vq^®Tl, (23) 

^vr = V^q®{q-q)+V^^®TX, (24) 

where g, q and vr are quark, antiquark and pion nonsinglet distributions , respectively, 
defined as 



q{,x) = qu{x) - qd{x) (25) 
q{x) = qu{x) - qd{x) (26) 
Ti{x]t) = 7r+(x;t) - 7r~(a;;t), (27) 

and, at one loop, the splitting functions and "P^^ are obtained from linear com- 
binations of the diagrams of Fig. 3 corresponding to various pseudoscalar mesons. 

Combining Eq.s (P^)-(P^ it is easy to show that that the evolution equation for 
quark distribution becomes 

(q + q)® {Qqq + Qqq) , (28) 
{q-q)® {Qqq - Qqq) + ® 27r, (29) 

where the various distributions and splitting functions are those which appear in 
Eq.s (p2D-(p7D. It follows that the evolution equation of Gottfried sum Eq.(^) is 
formally unchanged, provided of course the Qqq and Qqq are the new splitting func- 
tions which include meson emission. Notice that the meson distributions and splitting 



d 

dt^' 



q) 
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functions tt, V-^q and VqT, do not directly contribute to the evolution of g + g, and thus 
of So, whereas they do contribute to the evolution of g — In order to avoid double 
counting the contributions to the QCD evolution determined in this guise should not 
be added on top of the usual perturbative evolution. Rather, at large enough 
(presumably not larger than ~ 10 GeV^) the nonperturbative evolution should 
flatten out, so that at some large scale the perturbative behavior is reproduced. 

The next required step is the computation of the splitting functions. Now, the 
whole picture of evolution as being described by a probabilistic process relies on 
the dominance of certain diagrams (ladder diagrams) which contribute to the deep- 
inelastic scattering process, and are then effectively resummed by solving the evolu- 
tion equations. This corresponds physically to the dominance of processes where the 
emitted particle (a gluon in the standard evolution equations, a meson in the present 
case) is quasi-collinear to the original one (Weizsacker- Williams approximation) and 
has been proven [pO|] to hold for gauge theories (both abelian and nonabelian) in spe- 
cific gauges in the Bjorken limit. These proofs can actually be also extendedpO|] to 
the case of theories with "pseudoscalar glue", i.e., to the case of coupling of fermions 
to point-like pseudoscalars. It is much less clear that such a picture should hold in 
the case of extended pseudoscalars: as a matter of fact, coupling to pseudoscalar 
glue for all values of would lead to scaling laws for the moments of structure 
functions which are not in agreement with the experimental data and the effects 
discussed here can only be important in a limited range, not so large that the 
quark-antiquark structure of the mesons is seen. 

Hence, we proceed to the computation of the splitting function by assuming this 
approximation to hold. We verify then that this assumption is consistent, in that 
either significant evolution occurs and the meson emission cross section is dom- 
inated by a kinematical range where the meson is effectively pointlike, or, if the 
non-pointlike region is providing an important contribution to the cross-section, the 
evolution is negligible. This must happen at a value of low enough that the know 
asymptotic saling properties are unaffected. The splitting functions are then given 
by 

n,..(^;^)]n = ^^."?(^;^)' (30) 

in terms of the total cross section a'^'^{x; t) for the process of Fig. 3, integrated over 
all transverse momenta and expressed in terms of the usual Bjorken variables. Using 
this definition in the evolution equation (pOD the evolution of Sq is immediately found 



^This means that the so-caUed SuUivan process, deep- inelastic scattering off quarks inside pions 
in the nucleon wave function, cannot contribute directly to measurements of Sg- This is in keeping 
with the naive observation that a pion, being made of an equal number of u and d quarks (if quarks 
and antiquarks are counted with the same sign) cannot contribute to Sg- Pion emission may, 
nevertheless, contribute indirectly to Sg by modifying the flavor structure of the sea. Thus whereas 
pion counting does not tell us anything about Sg (contrary to occasional naive statements) pion 
emission mechanism (sometimes suggested within effective models jl^) may lead to a contribution 
to Sg- 
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to be given by the multiplicative law 



SG{t)=A{t,to)SG{to), (31) 
where the nonperturbative determination of A(t,to) is 



[^l(^'^o)]nonpert = ^Xp 



fT^' it) + ^af^' (t) - ^ant) 



(32) 



The problem of determining the evolution of Sq is thus reduced to the compu- 
tation of the various cross-sections which enter Eq.(^). The computation, which is 



somewhat involved technically, has been performed and discussed in Ref . |]18| . Be- 
cause the pion-quark coupling is of course unknown, the coupling has been written 
in the most general form compatible with Lorentz invariance; this requires introduc- 
ing four independent form factors (a pseudoscalar, two axial, and a tensor). The 
cross-section is then computed from the two diagrams (corresponding to s- and t- 
channel contributions) of Fig. 4. Using the asymptotic behavior of the form factors 
for large values of their arguments which is known from the behavior of the pion 
vertex function, it is found that for intermediate values of (up to a few GeV^) the 
cross sections are controlled essentially by the pseudoscalar form factor, while the 
strength of one of the axial couplings controls the large tail. The pseudoscalar 
form factor can be further constrained by use of the Ward identity and a condition 
which relates it to the physical value of the pion decay constant f^^. It turns out that 
all the evolution up to intermediate is thus parametrized by a single parameter 
which can be identified with the constituent quark mass, while an extra parameter 
(the axial coupling) controls the large tail. 

The main qualitative features of the result of this computation can be summarized 
as follows: 

1. The cross sections are indeed dominated by the region where the form fac- 
tor is pointlike; this is the region where the emitted meson and quark are 
almost coUinear, consistently with the parton interpretation. This is shown in 
Fig.5, where the cross section for pion emission is separated into its various 
pieces (Fig. 5a), which are then displayed (Fig. 5b, 5c) for various values of 
as a function of x. The pointlike dominance is apparent by noting that the 
form factor in the s and t channel are, respectively, ifg = (f{^{2s — M^)), and 
= (^(i(M2 - 2i)) with 



nid Is? + rn? 



^^(P^) = V^T^; (33) 



and s = as a; — 1, while t = —T^mi in the collinear /cnern limit, 

implying that the collinear, pointlike limit is x — 1 in the s channel and x ^ 
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in the t channel. It is interesting to observe that Fig.s 5b and 5c show that the 
growth of the cross section with (which determines, according to Eg. (|30|) , 
the anomalous dimensions) is due to the growth of the collinear peak, which 
dominates the cross section for all larger than ~ 0.1 GeV^. In the very 
small region (Q^ < 0.1 GeV^) where the non-pointlike contributions are 
important the dependence of the cross-section flattens (Fig. 5a). 

2. The axial coupling controls only the high-Q^ tail of the cross-section (see 
Fig. 5a), while the bulk of the dependence of the cross section occurs in a 
region (Q^ < 5 GeV^) where the cross section is controlled by the pseudoscalar 
coupling. 

3. The anomalous dimensions are positive definite throughout the range. The 
anomalous dimension for pion emission is indeed much larger than that for 
eta emission, thus leading to screening of the Gottfried sum, according to 
Eq.(0),(^). This is displayed in Fig.s 6a and 6b where the anomalous di- 
mension calculated from Eg. (pO|) is shown in the extreme cases of vr'' and t]' 
emission. 

4. All anomalous dimensions flatten at low (Fig. 6), thus so does necessarily 
the evolution of 5*^. It follows that this can be smoothly connected at the low 
scale where flattening occurs to a quark model predictions. 

5. All anomalous dimensions flatten at large (Fig. 6), due to the fact that the 
cross section for meson emission behaves as an inverse power of when is 
large enough. The evolution of 5*^ flattens more rapidly than any of the various 
evolutions due to meson emission separately does, since it is proportional to 
a difference of mesonic cross sections. Eg . (|3^) . It follows that the mesonic 
effects disappear at large Q^, where they behave as higher-twist effects, and 
perturbative behavior is regained. 

We proceed thus to a computation of the evolution of Sq according to Eq.s (pl|)- 
(P^. The starting value of Sg is chosen as the naive quark model expectation Eq.(|TD|). 
The scale at which this value holds is fixed either as \Q\ = 200 MeV, or as the scale 
at which the evolution flattens. Interestingly, this leads to the same result (small 
variations of the flattening point being irrelevant due to the multiplicative nature 
of the evolution). The final result for the evolution of So is shown in Fig. 7, for 
different values of the constituent quark mass (Fig. 7a) and for different values of 
the axial coupling which controls the high-Q^ tail (fig. 7b). This shows that the 
experimental value is reproduced as a consequence of nonperturbative evolution from 
the naive quark model value. In the region where the data are taken or slightly 
above the nonperturbative evolution flattens away; our computation provides an 
effective smooth interpolation (which is however sensitive to a free parameter of the 
computation, namely, the axial coupling) between the nonperturbative evolution and 
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the perturbative behavior which holds in the high hmit. A robust prediction of 
this computation is that significant evolution (characterized by anomalous dimensions 
of several orders of magnitude larger than the perturbative one) is taking place in 
the intermediate range, between roughly 0.1 and 10 GeV^. In particular, for 
reasonable values of the axial coupling one predicts that the asymptotic value of Sq 
is yet smaller than the presently reported one. 

In conclusion, nonperturbative evolution of the Gottfried sum due to bound-state 
emission in an intermediate energy range 0.3 GeV^ < \Q\ < 3 GeV^ is capable of re- 
producing the difference between the quark model value and the experimental result 
taken at = 4GeV^. The quark model value holds at a low scale, and is reduced by 
nonperturbative generation of a flavor-asymmetric sea, generated through a nonper- 
turbative generalization of the Altarelli-Parisi equations. The required asymmetry 
follows from the low-energy symmetry structure of QCD, as reflected in the spectrum 
of pseudoscalar mesons, and is due to the anomalous breaking of flavor singlet chiral 
symmetry. In this sense, the data on the Gottfried sum are displaying an effect of the 
axial anomaly. Because the nonperturbative evolution flattens, the value of Sq thus 
obtained is preserved by further QCD evolution: thus, nonperturbative QCD effects 
originating in the intermediate-energy region leave a spur even in the asymptotic 
limit. 



3 The sigma term puzzle 

We are flnally left with assumption 3, Eq.(^). The current knowledge of the nucleon's 
structure functions is in good agreement with this expectation: for example, at a 
scale of Q^=4 GeV^ strange quarks carry about 2-3% of the nucleon's momentum [21| 
(and accordingly much less at the nucleon's scale, since singlet anomalous dimensions 
are generally large and in this case threshold effects should also be relevant). This 
suggests that the matrix element {p\ss\p) should be roughlyf] of order 1 % of {p\uu\p). 
However, there is an old piece of data which seems to contradict this conclusion; this 
is known as the puzzle of the pion-nucleon sigma term0. The pion-nucleon sigma 
term is deflned as the nucleon matrix element of the light quark mass term in the 
QCD Hamiltonian: 

a = {N\Th{uu + dd)\N) (34) 

where rh is the average light quark mass, i.e., rh = {rriu + md)/2. The value of this 
operator can be extracted from experiment in two different ways; it would seem that 
the only way the results of the two determination can be made consistent with each 
other is by violating assumption 3. We show now, however, that a more careful 
analysis [^] reveals that this is not necessarily the case. 



^ It can actually be shown rigorously by defining parton distributions in terms of QCD light- 
cone operators that the forward matrix element of the operator ipiipi is to leading perturbative order 
equal to the first moment of the distribution of quarks plus antiquarks of flavor i. 
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The most direct way (at least theoretically) of determining the value of the sigma 
term is to observe that scalar quark densities appearing on the r.h.s. of Eq.(|3^) 
can be obtained by commuting a pseudoscalar charge with a pseudoscalar density 
according to the basic current algebra relation [p3| 

(35) 



where ip are quark isospinors and r° are Pauli isospin matrices. The pseudoscalar 
current is in turn related to the pion interpolating field byp^ 



(36) 



where a labels the pion according to its isospin and f^^a and m^a are the pion's 
decay constant and mass. Using Eq.s (|35| )(|36|) the matrix element on the r.h.s. of 
Eq.(|3^) can be related through straightforward current algebra manipulations [| 
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to a pion-nucleon matrix element which, in turn, is related to the isospin-even pion- 
nucleon scattering amplitude at vanishing momentum transfer. This, up to several 
subtle details of the extraction of the matrix element from the data, provides a direct 
unambiguous determination of a as defined by Eq. (|3^) , with the result EI 



a = 45 ±10 MeV. 



(37) 



On the other hand, if we assume Eq.(|) (i.e. assumption 3) to hold, we can equate 
the singlet and SU(3)-octet sigma terms: 



cr 



cTg = m{N\{uu + dd - 2ss)\N) . 



(3J 



The reason why such an identification is interesting is that erg is proportional to the 
matrix element of the octet portion Tig of the QCD Hamiltonian, which coincides 
with the octet portion of the mass term since SU(3) flavor symmetry is only broken 
by masses in QCD: 



0-8 



■{NlHslN) 



1 — rris/m 
Hs = - {ffi — nis) {uu + dd — 2ss 



(39) 
(40) 



Because in turn the matrix element of the octet part of the Hamiltonian provides 
SU(3) breaking to first order in perturbation theory, it can be measured from SU(3) 
mass splittings, for example in the baryon octet. 

Indeed, it is well known that in the nonrelativistic quark model one assumes the 
mass splittings in the baryon octet to be given by first order in perturbation by an 
operator Og which transforms as an octet under flavor SU(3) (identified with the 
octet component of the strong Hamiltonian). Group theory then leads to 124 



Mb = Mo + Mi{B\Y\B) + M2{B\ 



/(/ + l)--r 



\B), 



(41) 
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where I and Y , respectively, are isospin and hypercharge operators. Eq.(^Tj) expresses 
the masses of the all the octet baryons in terms of the three free parameters Mj and 
is in excellent agreement with experiment; in particular the Gell-Mann-Okubo mass 
formula, which is a consequence of it, is one of the great successes of the quark 
model. The nucleon matrix element of the symmetry breaking operator can also 
be determined from Eq.(^TD using SU(3) algebra: 

{N\0^\N) = Ml- M2/2 = Ma - Mh = (2Mjv - Ms - Ms) /3. (42) 

Identifying thus the symmetry breaking operator Os with the octet component of the 
QCD Hamiltonian Tig, and using the known value of the (current) quark mass 
ratio nis/fh = 25 ± 5, Eq.s (p9D,(^2D determine the octet sigma term as erg = 25 ± 5 
MeV, in blatant contradiction with the value of a Eq.(^), i.e. with the identification 
of a and erg Eq.(^). This discrepancy is known as the sigma term puzzle. 

The most simple-minded conclusion, namely that erg is much smaller than a just 
because the nucleon matrix element of ss is large leads to rather dramatic (and hard 
to believe) consequences: because the strange quark is so much heavier than the light 
ones, it would follow that about one third of the nucleon's mass is carried by strange 
quarks, in contradiction to the parton model results discussed above. Alternative 
explanations, based on a failure of first-order perturbation theory (in which case it's 
Eq.(^) and its consequences which do not hold) are extremely hard to reconcile with 
the success of the corresponding quark model formulas. It is interesting to observe 
that, once again, effective models of the nucleon have no difficulty in accommodating 



this state of affairs: in the Skyrme model ||2^ indeed assumption 3 is violated and 
CTg and a differ because of a large strange contribution, whereas in the bag model p8 
assumption 3 is valid but first-order perturbation theory fails. 



Upon closer examinat ion p2| , however, the argument given above is seen to con- 



tain a loophole: even though it is a true empirical fact that octet mass splittings are 
accurately described by the first order perturbative formula Eq.(^), how do we know 
that the operator Og which appears in this formula is the same as the operator Tig 



which is related to the sigma term by Eq. (|39D ? Of course, classically this is necessarily 



the case: mass splittings are given by the symmetry-breaking portion of the Hamilto- 
nian, and that's what Tig is. However in a second-quantized field theory this need not 
necessarily be the case: quantum corrections may change the symmetry structure, in 
general, hence also the pattern of symmetry breaking (through a quantum anomaly). 
Now, it is well known that this is precisely what happens for masses in the strong 
interactions: the masses which enter quark model formulas are constituent quark 
masses, which for light quarks differ by two orders of magnitude from the current 
quark masses which appear in the fundamental Lagrangian (and Hamiltonian) of the 
theory. Hence, what Og should be identified with is the octet piece of the constituent 
quark mass term, i.e., with the octet piece of the effective low energy Hamiltonian of 



QCD|29|, which is generally quite different from the fundamental QCD Hamiltonian. 



Thus we have the possibility of escaping the above dilemma if it just so happens that 
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the octet portion of constituent quark masses and current quark masses are different: 
in sucli a case Eq.(|38D still holds, along with first-order perturbation theory and its 
quark model consequences, but, at least in principle, octet splittings calculated from 
Og could be smaller than those calculated from Tis, which would explain the above 
discrepancy without invoking a violation of assumption 3. 

This explanation is actually less contrived than it may seem at first. We can see 
this by noting that actually the relationship between the matrix elements of and 
Tis is an exact consequence of the so-called conformal anomaly equation. This is 
an expression for the trace of the energy-momentum tensor T'^'^ which is true at the 
operator level in QCD: 

T% = (1 + 7r.) E ^^^^^^ + ^'^rc;:., (43) 

where the sum runs over all quark flavors, G^'^ is the gluon field strength, P = ^jf^ 
is the beta- function for the strong coupling = and 7m (/w) = is the 

mass anomalous dimension. The last term on the r.h.s. of Eq.(^3|) is due to the 
conformal anomaly; the anomalous dimension 7^ is present because the anomaly 
term and the mass term in Eq. (|43|) are not separately scale invariant (i.e., they mix 
upon renormalization), while the energy momentum tensor is (up to surface terms). 
Classically, both 7^, and the gluonic contributions would be absent. 



Taking the baryon matrix element of Eq.(^3D and using Lorentz invariance [23 
(i.e., essentially the fact that the trace of the energy- momentum tensor is a measure 
of the violation of scale invariance) it follows that the baryon mass Mb is given by 

{B\ (1 + ^rn)Y.^^i^^^^ + ^GTG^Ji?) = M^. (44) 

This shows explicitly that indeed mass splittings are linear in the expectation value of 
a certain operator, as in lowest order perturbation theory: because of the conformal 
anomaly, the linear "approximation" is actually exact; quantum corrections appear 
as an additive contribution to the naive first order perturbative formula. However, 
the operator which ought to be used is that on the r.h.s. of Eq.(|^). Otherwise 
stated, Eq.(|l) shows that 



(45) 

(8) 



where the subscript (8) indicates the SU(3) octet component of the operator in square 
brackets. 

It follows that the baryon mass Mb is in general given by the sum of a current 
contribution Mb and a dynamical contribution Mg '■ 

Mb = M^ + ME; (46) 
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= {B\Y^m.^.^.^B), (47) 

i 

ME = (i?i^GrG;:.+7™E^^^^^^i^)- (48) 

^ i 

The octet portion of Tis provides the current mass sphttings of M^, while the octet 
portion of Og provides the full splittings of Mb and it reduces to Tts only in the limit 
in which both 7^ and the gluonic term can be neglected. It should be stressed that 
this statement is an exact consequence of QCD. 

Superficially, one may think that the gluonic operator in Eq. (|45| ) should be flavor 
singlet and thus should not contribute to mass splittings. In this case, one would 
be left with 7^, which is however not quite negligible, since it may be computed in 
perturbation theory with the result [p^ jm[Q = 1 GeV] = 0.27.| This, however, is 



not possible since, as already mentioned, the gluonic contribution is not separately 
renormalization-group invariant, thus, by mixing with the fermionic component it 
acquires a nonsinglet portion. Indeed, Eq. (^6D is the unique renormalization-group 
invariant decomposition of the baryon mass. 

In fact, we can give two separate qualitative reasons which suggest instead that in- 
deed should be anticorrelated to M^. The first reason is based on the observation 
that a Ward identity argument implies that the one-meson reducible contributions to 
and must cancel: 

(Mg)"™" + {MET""' = 0. (49) 



Here are defined as the contribution to the matrix elements Eq.(^)-(|15|) where 

the operators which enter the definitions Eq. (|47|) (|48|) couple to the baryon state by 
first coupling to a meson which then couples to the baryon. Because such contri- 
butions are expected (by pole-dominance arguments) to provide the bulk of , 
Eq. (|49|) suggests that has a large component anticorrelated to M^. 



An independent reason is given by the observation that Eq.(^4D can be viewed 



as a relation between current masses (given by Eg . (]47|) assuming that the ma- 
trix element of ipiipi is just the number of quarks plus antiquarks of flavor i) and 
constituent quark masses nig (which add up to the baryon mass Mb); accordingly, 
the decomposition Eq.(^6D determines the dynamical contribution to the constituent 
quark mass (i.e., the difference between constituent and current mass) as a function 
of the current mass mp: 



= mq{m) — = m'^{m). (50) 

A baryon mass splitting will coincide with the mass splitting computed at the current 
level, i.e., neglecting the dynamical contribution M^ Eq.(^S]), only if Eq.(|5(]|) 

^Notice that if this were the only effect of the conformal anomaly, the discrepancy discussed 
above would be even larger, since this extra term would make larger than T-^g, rather than 
smaller as required to explain the puzzle. 
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turns out to be a flat function, i.e., not to depend on tlie current mass, so tliat 
tlie current and constituent masses differ by a fixed, flavor independent amount. 
However, the quahtative behavior of the function m^{m) Eq. (|SD|) is actually known. 
Indeed, when the current mass vanishes, the constituent mass reduces to roughly 
300 MeV, i.e., a third of the mass of the nucleon, which is made of quasi-massless 
quarks; hence m^{m = 0) = 300 MeV. When the current mass tends to infinity, 
quarks become quasi-static, and the current and constituent masses coincide; hence 
limfn-,oo (fn) = 0. It follows that is a decreasing function of the constituent 
mass, and not a flat function. This implies that constituent mass splittings are 
in general smaller than current mass splittings (so, for instance, the strange-light 
splitting is smaller for constituent masses than it is for current ones): the contribution 
of MJ* to mass splittings is anticorrelated to that of Mj^. 

Both these qualitative arguments can actually be turned into quantitative (or at 
least semi- quantitative) ones[^]. Firstly, one can actually estimate (Mj')°™'' by pole- 
dominance arguments, thus obtaining a rough estimate of the value of M^. Also, 
it is possible to actually compute the function m^{m) by studying the full quark 
propagator: more specifically, it is possible to attempt an approximate resolution of 
the Schwinger- Dyson equation satisfied by the full quark propagator. This leads to 
a determination of the quark self-energy, thus also of the quark constituent mass; 
this in turn can be assumed to give the baryon mass by an additive relation (i.e., the 
baryon mass is given by the sum of the masses of its constituent quarks). In either 
case, it is possible to determine AM, defined as 

AM= (A^IOslA^) - (A^l^sl^V)- (51) 

Even though the former determination is rather crude, and the latter is affected by 
an underlying theoretical uncertainty due to lack of knowledge of the QCD dynamics 
in the infrared (as reflected, for example, by the behavior of the strong coupling in 
the infrared limit), both determinations lead to a value AM ~ 150 MeV (with large 
uncertainties, perhaps up to 50 %). Specifically, with@ 30 MeV < AM < 250 MeV, 
and using Eq.(pID to relate Eq.(^TD (mass sphttings) to Eq. ([39|) (matrix elements of 
the sigma term) leads to the value 30 MeV < (Tg < 60 MeV, in perfect agreement 
with the value of a Eq.(^), i.e. with the identification of a and erg, consistently with 
assumption 3 [Eq.(^)]. 

In conclusion, if one takes into account that even "gluonic" quantities such as the 
dynamical contribution to baryon masses can have a flavor nonsinglet component, it 
follows that mass splittings at the constituent and current level are not the same. 
This allows one to resolve the puzzle of the discrepancy between the two different 
determinations of the sigma term, without invoking a large strange content of the 
nucleon, and without violating the linear quark model mass formulas, such as those 
leading to the Gell-Mann-Okubo relation. 
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4 Conclusions 



There are two main lessons to be learnt from this work, a phenomenological one 
and a theoretical one. The phenomenological lesson is that instances of discrepancy 
between experimental data and expectations based on the quark model are resolved 
once nonperturbative effects are taken into account. The theoretical lesson is that 
these effects appear both in the QCD evolution of operators, and in the computation 
of the operators themselves (in particular, in the determination of their symmetry 
structure). Such contributions, even though of course they originate in the gluon 
dynamics, depend strongly on the quantum numbers of the physical states under 
investigation and cannot be simply parametrized by universal flavor- singlet param- 
eters, such as vacuum condensates. These two points, taken together, suggest that 
new data may actually lead us to a better understanding of the infrared dynamics of 
QCD and thus to go beyond the very successful but limited set of predictions which 
perturbative QCD allows to make. 

Acknowledgements: I thank the organizers of the school and especially 
M A. Nowak for managing to assemble an unusually stimulating set of talks 
and discussions. This paper is largely based on work done in collaboration with 



R. D. Ball|[T^[^; I also thank V. Barone for discussions. This paper is partially 
based on notes taken by M. Engelhardt and K. Lucke (University of Erlangen) at a 
lecture given by me at the Erlangen-Regensburg Graduiertenkolleg. 



17 



References 

[1] J. Ashman et al., Phys. Lett. B206 , 364 (1988); Nucl. Phys. B328, 1 (1990). 

[2] For a review see G. Altarelli, in "The Challenging Questions", Proc. of the 1989 
Erice School, A. Zichichi, ed. (Plenum, New York, 1990). 

[3] Phys. Lett. B224 , 189 (1989); S. Forte, Nucl. Phys. B 331, 1 (1990); S. Forte 
and E. V. Shuryak, Nucl. Phys. B 357, 153 (1991). 

[4] P. Amaudruz et al., Phys. Rev. Lett. 66, 2712 (1991). 

[5] P. Amaudruz et al., CERN preprint CERN-PPE/93-117 (1993). 

[6] T. Sloan, G. Smadja and R. Voss, Phys. Rep. 162 , 45 (1988). 

[7] For a review see e.g. R. L. Jaffe and C. L. Korpa, Comm. Nucl. Part. Phys., 17 
163 (1987). 

[8] G. Altarelli, Phys. Rep. 81, 1 (1982). 

[9] See e.g. J. C. Collins, D. E. Soper and G. Sterman, in "Perturbative Chromo- 
dynamics", A. H. Mueller ed. (World Scientific, Singapore, 1989). 

[10] P. Amaudruz et al, Nucl. Phys. B371, 3 (1992). 

[11] see e.g. E. Leader and E. Predazzi, "Gauge Theories and the 'New Physics' " 
(Cambridge University, Cambridge, 1982). 

[12] K. Gottfried, Phys. Rev. Lett. 18, 1174 (1967) 

[13] V. R. Zoller, Phys. Lett. B279, 145 (1992); B. Badelek and J. Kwiecihski, Nucl. 
Phys. B370, 178 (1991). 

[14] S. Forte, Phys. Rev. D47 , 1842 (1993). 

[15] M. Anselmino and S. Forte, Torino preprint DFTT 6/92 (1993), Z. Phys. C, in 
press. 

[16] D. A. Ross, C. T. Sachrajda, Nucl. Phys. B 149, 497 (1979). 

[17] E. J. Eichten, 1. Hinchliffe and C. Quigg, Phys. Rev. D45, 2269 (1992). 

[18] R. D. Ball, S.Forte, Oxford and Torino preprint OUTP-93-18P and DFTT 9/93 
(1993). 

[19] A. W. Thomas, Nucl. Phys. A532, 177 (1991); S. Kumano, Phys. Rev. D43, 59 
(1991); D43, 3067 (1991); S. Kumano and J. T. Londergan, Phys. Rev. D44,717 
(1991). 



18 



[20] See e.g. Y. L. Dokshitzer, D. I. Dyakonov and S. I. Troyan, Phys. Rep. 58, 269 
(1980) 

[21] See V. Barone et al, Zeit. Phys. C58, 541 (1993); A. D. Martin, W. J. Stirling 
and R. G. Roberts, Phys. Rev. D47, 867 (1993). 

[22] R. D. Ball S. Forte and J.Tigg, Oxford and Torino preprint OUTP-92-35P and 
DFTT 69/92 (1993) 

[23] See e.g. V. de Alfaro, S. Fubini, G. Furlan and G. Rossetti, "Gurrents in Hadron 
Physics" (North-Holland, Amsterdam, 1973). 

[24] T.-P. Ghcng and L.-F. Li, "Gauge Theories of Elementary Particle Physics" 
(Glarendon, Oxford, 1984). 

[25] J. Gasser, H. Leutwyler and M. E. Sainio, Phys. Lett. B253, 252 (1991). 

[26] See e.g. J. Gasser and H. Leutwyler, Phys. Rep., 87 77 (1982). 

[27] J. F. Donoghue and G. R. Nappi, Phys. Lett. 168B, 105 (1986). 

[28] R. L. Jaffe, Phys. Rev. D21, 3215 (1980). 

[29] See e.g. H. Georgi "Weak Interactions and Modern Particle Theory" (Benjamin- 
Gummings, Menlo Park, GA, 1984). 



19 



Figure Captions 

Fig. [1] The gluon radiation diagram which leads to one-loop Altarelli-Parisi evolu- 
tion. 

Fig. [2] The two-loop diagram which generates nonsinglet evolution of quark distri- 
bution. When i = j the final state must be antisymmetrized with respect to the two 
identical quarks. 

Fig. [3] The meson radiation diagram which generates nonperturbative evolution. 

Fig. [4] Deep-inelastic scattering off a quark which radiates a bound state 11: a) 
t-channel diagram; b) s-channel diagram. 

Fig. [5] The 7r° emission cross section (with constituent quark mass M = 325 MeV) 
computed from the two diagrams of Fig.4: a) cross section integrated over x as a 
function of Q^; full line: full cross section with no axial coupling; dot-dash line: 
full cross section with axial coupling; dashed line: t-channel contribution (no axial 
coupling), diagram Fig. 4a; dotted line: s-channel contribution (no axial coupling), 
diagram Fig. 4b. b) t-channel contribution (dashed line of Fig. 5 function of x 

and Q^. c) s-channel contribution (dotted line of Fig. 5a) as a function of x and Q^. 

Fig. [6] Anomalous dimension calculated from Eq.(pOD for meson emission, for dif- 
ferent values of the constituent quark mass: a) 7r° emission; b) t]' emission. The axial 
coupling is set to a large (i.e. maximal) value. 

Fig. [7] Scale dependence of Sq computed from Eq. (|3TD : a) variation with the con- 
stituent quark mass; the axial coupling is set to a maximal value; b) Mg = 350 MeV; 
full curve: maximal axial coupling; dashed curve: no axial coupling. The experimen- 
tal point Eq.([TT|) is also shown. 
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